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ABSTRACT 

Superbursts are hours-long X-ray flares attributed to the thermonuclear runaway burn¬ 
ing of carbon-rich material in the envelope of accreting neutron stars. By studying the 
details of the X-ray light curve, properties of carbon combustion can be determined. In 
particular, we show that the shape of the rise of the light curve is set by the the slope of 
the temperature profile left behind by the carbon flame. We analyse RXTE/PCA ob¬ 
servations of 4U 1636-536 and separate the direct neutron star emission from evolving 
photoionized reflection and persistent spectral components. This procedure results in 
the highest quality light curve ever produced for the superburst rise and peak, and 
interesting behaviour is found in the tail. The rising light curve between 100 and 1000 
seconds is inconsistent with the idea that the fuel burned locally and instantaneously 
everywhere, as assumed in some previous models. By fitting improved cooling models, 
we measure for the first time the radial temperature profile of the superbursting layer. 
We find d\nT/d\n P Ri 1 / 4 . Furthermore, 20% of the fuel may be left unburned. This 
gives a new constraint on models of carbon burning and propagation in superbursts. 

Key words: accretion, accretion disks — stars: neutron — stars: individual: 4U 
1636-536 — X-rays: binaries — X-rays: bursts 


1 INTRODUCTION 

In the past decade rare long (hours-day) X-ray flares have 
been observed from neutron stars that accrete from a com¬ 
panion star in low-mass X-ray binaries (Cornelisse et al. 
2000). In contrast to the shorter (10-100 s) Type I X- 
ray bursts powered by hydrogen/helium flashes (e.g., 
Lewin et al. 1993), these so-called superbursts are thought 
to arise from runaway carbon fusion in the ashes of hy¬ 
drogen and helium burning (Strohmayer & Brown 2002; 
Gumming & Bildsten 2001). 

Gumming & Macbeth (2004) assumed that after super¬ 
burst ignition the carbon would rapidly undergo complete 
burning to iron group elements, and modelled the super¬ 
burst light curve by numerically following the subsequent 
cooling of the neutron star envelope. These cooling models 
were fit to superburst light curves to derive the depth and 
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energy-content of the burning layer (Gumming et al. 2006). 
Although they are a good match to the tail of superbursts, 
these models failed to reproduce the rise of the observed light 
curves. In particular, the cooling model light curves predict 
that the luminosity decreases between 100 to 1000 seconds 
after onset of the superburst, whereas the observed light 
curves are increasing in luminosity at that time. With the 
same assumption of complete burning, Weinberg & Bildsten 
(2007) also found declining light curves at 100-1000 seconds. 
In contrast, Keek & Heger (2011) carried out fully self- 
consistent simulations of superbursts that followed the car¬ 
bon burning and energy deposition in detail. They instead 
found a rising luminosity during this phase of the superburst 
in some models (where the precursor did not dominate the 
rise). 

This earlier work indicates that the phase of the super¬ 
burst light curve prior to the peak at 1000 seconds is sens¬ 
itive to the details of how the carbon combustion deposits 
energy in the neutron star envelope. The physics of how the 
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carbon burning propagates is complex and may not be accnr- 
ately captured by current ID models. Rather than explicitly 
treating the hydrodynamic process of flame spreading, these 
models resolve the burning layers in the radial direction and 
include a ID approximation to turbulence. In this paper we 
derive a useful constraint on such modelling by using the ob¬ 
served light curve to constrain the radial temperature profile 
left behind by carbon burning. This may inform other ap¬ 
plications of carbon flames, such as the study of detonation 
and deflagration in Type la supernovae (e.g., Woosley et al. 
2011 ). 

Most of the current 24 (candidate) superbursts (e.g.. 
Keek et al. 2012; see Negoro et al. 2012 and Serino et al. 
2014 for recent detections) are detected with instruments 
that have limited sensitivity (e.g., in’t Zand et al. 2003; 
Keek et al. 2008). Furthermore, the rise is often not ob¬ 
served due to gaps in the data from Earth occulta- 
tion or passage through the South-Atlantic Anomaly (e.g., 
in’t Zand et al. 2004). The two highest quality observations 
were performed with the Proportional Counter Array (PCA; 
Jahoda et al. 2006) on-board the Rossi X-ray Timing Ex¬ 
plorer {RXTE-, Bradt et al. 1993). The first was a bright 
superburst from 4U 1820-30 (Strohmayer & Brown 2002; 
Ballantyne & Strohmayer 2004). Its luminosity quickly 
reached the Eddington limit, which impedes the study of 
the intrinsic rise (e.g., in’t Zand & Weinberg 2010). The 
second PCA superburst observation was from 4U 1636-536 
(Strohmayer & Markwardt 2002; Kuulkers et al. 2004). It 
was not Eddington-limited, and at present it provides the 
only detailed observation of a superburst rise. Moreover, it 
is the only (super) burst for which spectral analyses were 
able to separate the direct burst emission from an evolving 
persistent component as well as photoionized reflection off 
the accretion disc (Keek et al. 2014a, b). Reflection refers to 
scattering of burst emission off the disc, which introduces 
features in the spectrum, such as a strong Fe Kq emis¬ 
sion line, that depend on the ionization state of the disc 
(Ballantyne 2004). 

In this paper we fit improved cooling models to the su¬ 
perburst light curve of 4U 1636-536. First we investigate the 
robustness of the light curve derived from the observations, 
in particular with respect to the influence of deviations of the 
neutron star’s spectrum from a blackbody (Section 2). Sub¬ 
sequently, we create improved cooling models to constrain 
properties of the superbursting layer, including the slope of 
its temperature profile (Section 3). Comparison of the super¬ 
burst tail to the best-htting model light curve reveals inter¬ 
esting, but not entirely understood, behaviour (Section 4). 
Finally, we discuss how carbon combustion produces the in¬ 
ferred temperature gradient (Section 5). 


2 RECONSTRUCTING THE OBSERVED 
LIGHT CURVE 

For comparison to super burst cooling models, we set out 
to determine the most reliable light curve of neutron-star 
emission by separating out emission from and reflection off 
the accretion disc. In previous analyses, a blackbody was 
used to model the thermal spectrum of the neutron star 
(Kuulkers et al. 2004; Keek et al. 2014a,b). A neutron star 
atmosphere spectrum is, however, expected to exhibit de¬ 


viations from a blackbody due to, e.g., electron scatter¬ 
ing which influence the blackbody parameters obtained in 
a spectral analysis (Suleimanov et al. 2011). To investigate 
how this affects the inferred net superburst light curve, we 
repeat the analysis procedure of Keek et al. (2014b) with 
neutron star atmosphere spectra instead of a blackbody. 

2.1 Observations and spectral models 

We reuse the PCA source spectra, instrumental back¬ 
ground spectra, and response matrices that were created for 
Keek et al. (2014b) (see also Keek et al. 2014a). The spec¬ 
tra are collected in 64 s time intervals, and we consider the 
3—20 keV energy range in our analysis. The superburst spans 
4 RXTE orbits. As the signal is weaker in the final orbits, 
we use a single spectrum covering orbit 3. Orbit 4 is omit¬ 
ted, because the reflection features cannot be detected, and, 
therefore, the reflection signal cannot be separated from dir¬ 
ect burst emission (Keek et al. 2014b). Spectral fitting is 
performed with XSPEC version 12.8.1 (Arnaud 1996), and 
1 a errors in the fit parameters are derived using a change 
in goodness of fit of Ax^ = 1. 

The spectra are described using 3 components; persist¬ 
ent emission from the accretion disc, direct thermal emission 
from the neutron star, and reflection of the thermal emission 
off the disc. The persistent part is modelled by a power law 
with an exponential high-energy cut-off at energy Ecutoft- 
For the other two components Keek et al. (2014b) employed 
a blackbody with temperature kT and a model of a reflec¬ 
ted blackbody (Ballantyne 2004) at the same temperature, 
as a function of the disc’s ionization parameter ^ (in units of 
erg s“^ cm in this paper). We refer to this model as BB. The 
spectral model further includes the photoelectric absorption 
model vphabs with an equivalent hydrogen column Nh (c.f., 
Pandel et al. 2008) and smoothing of the reflection compon¬ 
ent due to Doppler broadening and relativistic effects using 
the rdblur model (Fabian et al. 1989). 

As an alternative for the blackbody, we use neut¬ 
ron star atmosphere spectra based on the models by 
Suleimanov et al. (2012) with a surface gravity of logg = 
14.3 (g in units of cms“^). Instead of the blackbody tem¬ 
perature, these models are a function of the ratio of the 
luminosity to the Eddington limited luminosity: L/LEdd- 
We create two grids of atmosphere spectra each containing 
100 values of the Eddington ratio, equidistant on a logar¬ 
ithmic scale, in the range 10“® < L/LEdd < 1.08. These 
grids are optimized for the analysis of the superburst, as 
they have improved resolution at lower values of L/Lsad 
compared to Suleimanov et al. (2011). One grid is for solar 
composition (model SI), and the other has a 100 times lower 
metal content (model SOOl). During the spectral fits with 
each grid, we fit for L/LEdd and the normalization. For both 
grids we calculate photoionized reflection spectra using the 
code of Ballantyne et al. (2002) using the same procedure as 
Ballantyne (2004). 

2.2 Spectral analysis 

Because of the limited data quality. Keek et al. (2014b) 
found that certain parameters need to be constrained dur¬ 
ing the spectral hts. In particular, in the hrst RXTE orbit 
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Figure 1. For 3 spectral models, flux as a function of time, t, for the direct neutron star emission (F^g), reflection (Fjefl), and persistent 
flux (Fpers; the long-dashed line is the pre-superburst value). At the bottom we show the reflection fraction Fred/P)sis- Fns is modelled 
with a blackbody in BB and by a neutron star atmosphere in SI and SOOl. SI assumes solar metallicity, whereas SOOl uses a 100 times 
lower metalicity. The left and right panels each have a different horizontal scale. The Icr uncertainty is only indicated for the last data 
point. Excellent fits are obtained around the peak (left panel), where Fp^g for the 3 models agrees at all times within 17%. Issues in the 
tail appear as large unphysical variations of F^efl for SI and SOOl (right panel; see also Section 5.4). 


(excluding the first 128 s, which contain the pre-superburst 
flux and the precursor burst; see Figure 1) all spectra are 
fit simultaneously, and the normalizations of the power law 
and the direct neutron star component as well as A^h are 
tied between the spectra. The power-law index and Fcutoff 
are fixed to the pre-superburst values (Keek et al. 2014a). In 
orbits 2 and 3, Fcutoff is free, but the normalization of the 
neutron star component is fixed to the value obtained for the 
first orbit. This procedure yielded a good fit with model BB 
(Keek et al. 2014b), and repeating it with model SI gives a 
similar goodness of fit. It does not, however, result in a good 
fit for SOOl in the first orbit. To obtain a good fit with SOOl, 
we allow Fcutoff to vary in the first orbit, keeping its value 
tied between the spectra in the first orbit. The best-fitting 
value, Fcutofi = 6.17 ± 0.07, is significantly larger than the 
pre-superburst value of Fcutoff = 4.80 ± 0.06 (Keek et al. 
2014a). The xt values obtained with SI and SOOl have a 
similar distribution as the values for model BB (Keek et al. 
2014b). 

SI and SOOl yield parameter values that are over-all 
quantitatively consistent with model BB. In particular, the 
weighted mean of log^ in the first orbit is within 2.8% of 
the BB value of log^ = 3.384 ± 0.009, and in the second 
orbit all models exhibit a transition to a lower value of log^. 
Furthermore, the Eddington ratio of the atmosphere mod¬ 
els corresponds to an effective temperature of the neutron 
star photosphere that is lower than the measured blackbody 
temperature (Suleimanov et al. 2011). The ratio of the two 
is referred to as the colour correction factor, /c. For model 
SI, we find /c = 1.464 ± 0.007 when the Eddington ratio 


peaks at L/L-^dd = 0.499 ± 0.007, and it drops to a mean 
value of fc = 1.295 ± 0.010 in the first 512 s of the second 
orbit where L/L^dd = 0.101 ± 0.002. Similarly, for model 
SOOl the Eddington ratio peaks at L/L^dd = 0.441 ± 0.006 
with /c = 1.511 ± 0.008, and maintains a comparable value 
of fc = 1.483 ±0.011 while the Eddington ratio decreases to 
L/Lsdd = 0.0604T0.0013 at the start of the second orbit. In 
both cases the behaviour is consistent with the theoretical 
predictions for the considered values of L/L-Edd'- fc is close 
to constant for a metal-poor atmosphere, but it exhibits lar¬ 
ger variations for solar composition (Suleimanov et al. 2011, 
2012). The deviations of the spectrum from a blackbody are, 
therefore, too small for the data to discriminate between SI 
and SOOl, and the models simply behave as predicted. 

The neutron star fluxes, Fns, are consistent for BB, SI, 
and SOOl: in the first orbit within 17% (Figure 1). Because 
SOOl has a higher Fcutoff, its persistent flux, Fpers, is larger 
by 29% in the first orbit. In the following orbits SOOl’s Fpers 
is consistent with BB and SI. Furthermore, in the second or¬ 
bit the reflected neutron star flux, FreU, of both SI and SOOl 
exhibits strong variations. For some spectra the fit prefers 
low values of L/L^dd, which require large bolometric correc¬ 
tions and produce unrealistically large values for Feed and 
for the reflection fraction Frefi/FNS. For other spectra, the 
fits find more reasonable values of Frefl, and this ‘track’ is 
closer to the fluxes and reflection fractions obtained for BB 
(Figure 1). We attribute these issues to the increasingly poor 
data quality in the superburst tail. 

Because of the similarity of Fns for all models around 
the peak, and the issues of SI and SOOl in the tail, we choose 
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model BB for comparison with theoretical light curves. We 
calculate the net superburst luminosity from Fns assum¬ 
ing isotropic emission by the neutron star at a distance 
of 6.0 ± 0.5 kpc (Galloway et al. 2008). We do not include 
the distance error in the uncertainty of the luminosity data 
points, but it is included as a systematic uncertainty when 
we fit cooling models (Section 3). Because of the mentioned 
issues in the tail of the superburst, we first focus on the rise 
and the peak of the superburst in the first orbit. In Section 4 
we look at the tail in more detail. 


3 COOLING MODELS FOR THE 

SUPERBURST LIGHT CURVE 

In this section, we make models of the superburst light curve 
by assuming that the nuclear burning occurs rapidly after 
ignition, leaving behind a temperature profile that we then 
take as an initial condition for the subsequent cooling. We 
show that the light curve slope during the rising phase at 
100-100 s constrains the temperature prohle in the layer set 
by carbon burning. 

3.1 Details of the models 

To calculate the cooling of the heated layer, we follow the ap¬ 
proach of Cumming & Macbeth (2004) and Cumming et al. 
(2006). The main differences from those calculations are in 
the choice of initial temperature prohle for the cooling, a 
correction to the radiative opacities used in those papers, 
and a more accurate outer boundary condition. 

To follow the cooling of the layer, we integrate^ the heat 
equation 

dT dF ^ AacT^ dT 

where F is the heat hux, e,, the neutrino emissivity, T the 
temperature, y the column depth, k the opacity, a the radi¬ 
ation constant, and c the speed of light. We use the method 
of lines, in which we difference the spatial part of these 
equations on a grid in column depth with a constant spa¬ 
cing Aa: = Alny, and then integrate the stiff set of or¬ 
dinary differential equations that result forward in time us¬ 
ing an implicit integrator. The inner zone is set to be a 
factor of 1000 deeper in column than the ignition depth. 
We typically set the outer zone in our time-dependent cal¬ 
culation at a column depth yt — 10® g cm“^, which is a 
small fraction of the typical fuel column (< 10“®). We ad¬ 
opt as default neutron star parameters M = 1.4 Mq and 
R = 12 km giving a gravity gi 4 — cm s“® = 1.60 

and redshift 1 + z = {1 — 2GM /(?= 1.24, where 
g = {GM/R^){l + z). This choice of g is slightly smaller than 
what was used for the atmosphere spectra (Section 2.1), but 
there is no problem with consistency, because we will fit 
the cooling models to the light curve from spectral model 
BB. The calculation is Newtonian in that redshift variations 
across the thin layer are ignored. Instead, we set the grav¬ 
ity to be a constant and then redshift the time and flux to 
infinity to compare to observations. 

^ The code used in this paper is available at 
https://github.com/andrewcumming/burstcool. 


The microphysics input is as follows. For the equa¬ 
tion of state of the electrons, we use the fits of Paczynski 
(1983), except for the electron heat capacity which is a mod¬ 
ified version of the Paczynski (1983) formula that has the 
correct limits for non-degenerate and degenerate electrons 
(Schatz et al. 2003). The electron Fermi energy is calculated 
as needed according to Chabrier & Potekhin (1998). Cou¬ 
lomb corrections for the ions are calculated using the res¬ 
ults of Potekhin & Chabrier (2000) (see their Equation. 17). 
The free-free and electron scattering radiative opacities 
are calculated following Schatz et al. (1999). When form¬ 
ing the total opacity, we include an additional factor from 
Potekhin & Yakovlev (2001) to account for the fact that 
Rosseland mean opacities are not additive (see discussion 
in Stevens et al. 2014). The thermal conductivity is calcu¬ 
lated by implementing the htting formulas of Potekhin et al. 
(1999). Neutrino emissivities due to the plasma, pair, and 
bremsstrahlung processes are included using the fitting for¬ 
mulae from Schinder et al. (1987) and Haensel et al. (1996). 

The boundary condition F{T) at y = yt is determined 
from a grid of constant flux atmospheres which are integ¬ 
rated separately. The key region is the “sensitivity strip” at 
column depths ~ lO’^-lO® g cm“^, which is assumed to con¬ 
sist of ®®Ni in our default model (matching the composition 
of the cooling layer). The layers y < 10^ g cm“^ are pure 
helium (the composition at low density does not affect the 
F~T relation, being out of the sensitivity strip). 

3.2 The initial temperature profile 

The starting temperature profile for the cooling is set to be 
a power law with depth, Tf{y) = Tb{y/yb)°‘■ The normaliz¬ 
ation Th is chosen such that the average energy input is a 
specihed value Enuc = EislO'^® erg g“^, 

rvb fTf(y) 

/ dy cpdT = yUnuc. (2) 

Jyt Jt, 

Since Tf ^ Ti, the light curves are not very sensitive to the 
choice of U; for simplicity we set U = 2 x 10® K, constant 
with depth. 

The effect of the choice of the slope of the temperature 
profile a on the light curve is shown in Figure 2. The light 
curves shown have the same total energy release and depth, 
but different choices of a. We also show two other choices of 
initial temperature profile. One is to assume no heat trans¬ 
port, and burn the fuel locally and instantaneously. This was 
the assumption of Cumming & Macbeth (2004) and gives a 
profile with a « 1/8. The second is an adiabatic initial pro- 
Hle, with dluT/dlnP = Vad at each depth. This could res¬ 
ult if convection is able to efficiently mix the burning layer. 

Figure 2 shows that different choices of a result in dif¬ 
ferent light curve slopes at times 100-1000 s before the light 
curve peak. Smaller values of a < 0.21 lead to a luminosity 
that falls with time during this phase; larger values of a give 
a luminosity that increases with time. We see that a meas¬ 
urement of the slope of the cooling curve at early times is 
a direct measure of the slope of the temperature profile in 
the layer immediately following burning. This is explicitly 
shown in Figure 3 which shows the temperature slope a as 
a function of the slope of the light curve measured at 200 
seconds. 

The form of the relation between a and dInL/dInf 
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Figure 2. Model light curves (top panel) with the same total 
energy release and ignition depth, but different initial temper¬ 
ature profiles (bottom panel). The solid curves have T oc P“ 
with (top to bottom) a = 0.2, 0.225, 0.25 and 0.3. The dot- 
dashed curve assumes local deposition of energy (a si 1/8 as in 
Gumming & Macbeth 2004); the dashed curve has an adiabatic 
profile dlnT/d\n P = Vad . For all models the ignition column is 
2 X 10^^ gcm“^ and energy release Eis = 0.25. Crosses are the 
data points of the observation of 4U 1636-536, where time starts 
at the onset of the precursor (t = 72 s in Figure 1). 

shown in Figure 3 can be understood using the analytic es¬ 
timates for the early time light curve from Appendix A of 
Gumming et al. (2006). Assuming constant opacity for sim¬ 
plicity, they modelled the early time light curve by writing 
the flux F oc T*/y, where y is the depth that the cooling 
wave has reached at a time t. This depth is given by set¬ 
ting the thermal time-scale to the current time t = itherm oc 
y"^/^ jT'^, With a temperature profile of the form T oc j/“, we 
find dlnL/dlnf = —4(4a — 1)/(8 q — 7). This qualitatively 
agrees with Figure 3, but is not accurate in detail as we 
might expect from a simple constant opacity estimate. We 
find that d\n L / d\nt = 5(a — 0.21) is a good approximation 
to the numerical results, shown in Figure 3 as a dotted line. 

3.3 Comparison to the data 

The data for the 4U 1636-536 light curve is plotted in Fig¬ 
ure 2. We take the start time t = 0 at the onset of the 
precursor burst (e.g., Keek et al. 2014a), which corresponds 



Lightcurve slope d In L/d In t at 200 seconds 


Figure 3. The slope of the light curve at 200 seconds 
(d In L/d In t) as a function of the temperature slope in the layer at 
the onset of cooling T oc P“. The dotted line is the analytically- 
motivated fit dlnL/dlnt = 5(a — 0.21). 


to t = 72 s in Figure 1. As described by Gumming et al. 
(2006), the total energy release is constrained by the lumin¬ 
osity of the early part of the light curve (before the peak 
at ~ 10® seconds), whereas the depth of the layer is con¬ 
strained by the cooling time-scale associated with the tail 
of the burst (after 10® seconds). We find that an energy re¬ 
lease Eis ~ 0.25 and a depth y ~ 2 x 10®® g cm~^ give a 
good fit to the data. In addition, as can be seen in Figure 2, 
the data strongly constrain the temperature profile to have 
a = 0.25. The observed slope d\n L/d\nt is positive, ruling 
out the local burning temperature profile, but also not steep 
enough to match the adiabatic profile. 

The curves shown in Figure 2 are for a particular 
choice of neutron star mass and radius, composition of the 
burning products, and distance (Section 3.1). We investig¬ 
ate the robustness of our conclusions to these parameters. 
We use the emcee Markov Ghain Monte Garlo (MGMG) 
code (Foreman-Mackey et al. 2013) to explore the parameter 
space. In doing so, we fit the data for t < 1400 s only. This 
includes the peak of the light curve, but excludes most of 
the decay. As can be seen in Figure 2, the observed lumin¬ 
osity shows abrupt changes of slope during the decay phase 
which can bias the fits, since there are many data points in 
the decay part of the light curve. By including the data in 
the rise and peak only we are able to obtain good fits to the 
data (x2 ~ 1) and we still obtain a tight constraint on the 
allowed column depth. 

We find that the derived values of Eis, y and especially 
a are not very sensitive to the uncertainties in the other 
parameters. For a fixed distance, we find that changing the 
heavy element composition, for example from Fe to Si, or the 
mass and radius of the neutron star within typical values, 
e.g. 1.2 to 2.QMq and 10—13 km changes the derived energy 
and column depth by « 10 %, with much smaller changes in 
a. Ghanging the distance has a larger effect on the derived 
energy and column depth (Gumming et al. 2006). For ex¬ 
ample, with a prior on distance from 4 to 8kpc, we find 
a range of E\s ~ 0.24—0.32 and y « 1.5—4 x 10®® gcm“^. 
The distribution of allowed slopes is still rather narrow, with 
a = 0.242 ± 0.004. 
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4 THE TAIL OF THE SUPERBURST 

Although the cooling models are in excellent agreement with 
the rise and peak of the observed light curve, in the tail Ens 
decreases much faster than predicted (see the corresponding 
luminosity in Figure 2). As the normalization of the neutron 
star spectral component is kept fixed in the tail, the beha¬ 
viour of the flux is purely due to the temperature evolution. 
This implies exceptionally fast cooling of the neutron star 
atmosphere. 

As an alternative assumption for the normalization of 
the BB blackbody component in the tail, we fix the reflection 
fraction to the mean value in the first orbit: Erefl/ENS = 
0.665 ± 0.012 (model BBfrac). The physical implication of 
a constant reflection fraction is that the accretion geometry 
is unchanged during the superburst. The goodness of fit for 
BBfrac is similar as the other models. This produces a slower 
decay that closely follows the trend of the cooling model 
(Figure 4). To quantify how well BBfrac matches the cooling 
model, we take the ratio of Ens to the model flux. In orbit 
2, the weighted mean is 0.858 ± 0.009. This implies that the 
actual reflection fraction in the tail is 0.43±0.02. A decrease 
of the reflection fraction may be expected if the disc was 
puffed up by X-ray heating in the peak, and returned to a 
flatter geometry as the illumination was reduced in the tail 
(e.g., Ballantyne & Everett 2005; Lapidus & Sunyaev 1985). 

Model BBfrac has a similarly fast decrease of the black- 
body temperature as BB. The slower decline of the flux, 
therefore, implies that the temperature decrease is partially 
compensated by an increase in the normalization. In the final 
1024 s of the second orbit, the weighted mean of the ratio 
of Ens for BBfrac and BB is 5.2 ± 0 . 2 . As the blackbody 
normalization is proportional to the apparent emitting area, 
this corresponds to a radius increase with respect to the peak 
by a factor 2.27 ± 0.04. The appearance of radius expansion 
in this part of the superburst is puzzling, and in Section 5.4 
we further discuss the challenges of the interpretation of the 
superburst tail. 


5 DISCUSSION 

We first discuss the main result from our analysis: con¬ 
straints on the superburst parameters from cooling mod¬ 
els, including the temperature profile created by the carbon 
flame. Next we compare the inferred ignition column to the 
observational constraints on the accreted column. Finally, 
we discuss the robustness of the observed light curve and in 
particular the issues with the interpretation of the spectra 
in the superburst tail. 

5.1 Retracing the carbon flame with cooling 
models 

Carbon combustion has been studied extensively in the con¬ 
text of Type la supernovae, in particular to investigate 
whether carbon flames propagate as a detonation and/or a 
deflagration (e.g., Woosley et al. 2011). Superbursts provide 
a unique opportunity to study carbon flames, as they occur 
close to the surface of the star, without the thermonuclear 
event disrupting it. Weinberg & Bildsten (2007) calculated 
the propagation of a detonation through the carbon layer. At 


depths y < yaet ~ 2 x 10 ^^ g cm“^ the time-scale to burn 
carbon is longer than a dynamical time-scale, so that the 
detonation ceases and the shock runs ahead of the burning. 
The column depth that we infer for the 2001 4U 1636-536 
superburst is about this depth, suggesting that a detonation 
did not likely occur in this superburst. 

Our fits with cooling models find a temperature pro¬ 
file with power-law slope a ~ 1/4, whereas complete local 
burning produces a = i/s. One way in which a temperat¬ 
ure profile with a > i/s could come about is if the burn¬ 
ing becomes more and more incomplete at low densities. 
Weinberg & Bildsten (2007) pointed out that the shallow 
superburst light curves imply a steeper temperature gradi¬ 
ent than would result from a complete burn. They suggested 
that perhaps the convective deflagration leaves behind an in¬ 
creasing amount of unburned fuel as it moves to lower dens¬ 
ities. Since local instantaneous burning produces a temper¬ 
ature Tf oc (Cumming & Macbeth 2004), E/ oc E“ 

implies Eis oc p(8«-i)/4, p^j. q, = 1/4 this is Eig oc E^/^. If 
the fuel layer initially has a homogeneous composition, this 
describes the fraction of the fuel that is burned as a function 
of depth. For example, assuming that at the ignition depth 
100% was burned, we estimate that only 50% was burned 
at an order of magnitude shallower depth. Integrating over 
all depths, as much as 20 % by mass of the fuel may be left 
unburned. 

Sound waves or a shock generated by carbon igni¬ 
tion have been suggested to heat the atmosphere above 
the superburst fuel layer and induce hydrogen and helium 
burning, which would power the precursor bursts observed 
at the onset of superbursts (Weinberg & Bildsten 2007; 
Keek & Heger 2011; Keek 2012). The tail of these precurs¬ 
ors and subsequent stable burning of freshly accreted hydro¬ 
gen and helium may lower the steepness of the light curve 
(Keek et al. 2012). We have reduced this effect by starting 
our fits well after the precursor. 

Keek & Heger (2011) and Keek et al. (2012) calculated 
fully self-consistent Id multi-zone models of superbursts. 
The light curves show both rising and declining phases 
at times hundreds to thousands of seconds, depending on 
the accretion rates and base fluxes chosen for the models. 
For example, the superburst light curve shown in Figure 3 
of Keek et al. (2012) shows a rise between « 30 and 300 
seconds, and then a declining light curve from 300 to 10^ 
seconds. Our new cooling models are able to reproduce the 
light curves of these models. In a future publication, we will 
study in detail how the measured temperature profile is re¬ 
produced in these Id models. 

Given that carbon burning to nickel produces 
1.0 MeV u“^ = 0.96 X lO*^® ergg“^ (Schatz et al. 2003, ignor¬ 
ing a possible contribution from photodisintegration), the 
measured Eis translates to a carbon mass fraction of 26%. 
This is consistent with the value previously measured with 
local burn models (Cumming et al. 2006). Furthermore, we 
find an ignition column that is less than half as large as 
previously measured. This is due to the different shape of 
our new model light curves. Therefore, the inferred column 
depths for other superburst observations may have been 
too large as well (Cumming et al. 2006; Keek et al. 2008; 
Kuulkers et al. 2010; Altamirano et al. 2012). For fits to the 
light curve of the candidate superburst from EXO 1745- 
248, Altamirano et al. ( 2012 ) used an early version of our 
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Figure 4. Comparison of two derived superburst light curves to the cooling model that fits best to the first 1400 s (top of the shaded 
region; see Section 3). BB assumes a constant normalization of the spectral component that describes the thermal emission of the neutron 
star, whereas spectral model BBf^ac assumes the reflection fraction to be constant. The left and right panels and their scale are the 
same as in Figure 1. 


cooling models with an a value close to the one that we 
obtain, in anticipation of this paper. Because the onset of 
that superburst was not observed, it did not itself provide 
a constraint on a. The 2001 superburst from 4U 1636-536 
is the highest quality observation of the rise and peak, and 
improved measurements of superburst parameters including 
a, therefore, require new superburst observations with fu¬ 
ture large-area X-ray telescopes such as NASA’s Neutron 
Star Interior Composition Explorer {NICER, scheduled for 
launch in 2016; Gendreau et al. 2012). 

5.2 Accretion history 

Three other superbursts have been observed from 4U 1636- 
536. The one immediately preceding the super burst dis¬ 
cussed in this paper was 1.75 yr earlier (Kuulkers 2009; 
no superburst has been observed from this source after¬ 
wards). From the Multi INstrument Burst ARchive (MIN- 
BAR; Keek et al. 2010), which is the largest catalogue of 
X-ray observations of bursting sources, we select all 132 
RXTE PGA (Galloway et al. 2008) and BeppoSAX WFC 
(Cornelisse et al. 2003) observations in that time interval. 
The observations are spread out in time, and constitute a 
total exposure time of 1.44 Ms. The unabsorbed persistent 
3 — 25 keV flux varies between observations by a factor 2.8. 
We determine its time-averaged value as the mean of all 
flux measurements weighted by the duration of each obser¬ 
vation: Tpers = (4.17 ±0.07) X 10~® ergs“^ cm“^. Employ¬ 
ing a bolometric correction of 1.1 (Fiocchi et al. 2006) and 
a distance of 6.0 kpc, the average persistent luminosity is 
Tpers = (1.98T0.03) X 10®^ ergs“^ (distance uncertainty not 
included). For 1.75 years of accretion on to a neutron star 
with parameters as described in Section 3.1, this corresponds 
to an accretion column of j/acc = (3.53±0.05) x 10^^ gcm“^, 
assuming the efficiency of releasing gravitational potential 
energy in X-rays is 100%. Because of our assumptions, 
the actual uncertainty is probably at least several tens of 
percent. Given the best-fitting value for i/ign from cooling 
models (Section 3.3), j/acc can accommodate 2 superbursts: 
2 /acc/t/ign = 1.8±0.2, where we use the MGMC-derived error 
on j/ign without the uncertainty in the distance, which drops 
out. We regard it likely that an extra superburst occurred 


roughly in the middle of the preceding 1.75 years. Around 
this time no superburst has been detected with, e.g., the 
All-Sky Monitor on RXTE (e.g., Kuulkers 2009) or other in¬ 
struments, but this could easily be attributed to the sparse 
coverage of the source by the X-ray observatories that were 
operational at the time. 

5.3 Robustness of the observed rise and peak 

In the rise and around the peak of the superburst light curve, 
the spectrum clearly separates in three components, one of 
which we interpret as the direct thermal emission from the 
neutron star. The best-fitting parameters of the considered 
spectral models for this component can be somewhat differ¬ 
ent: e.g., the peak value of L/L^dd for the metal-poor at¬ 
mosphere model is 12% lower than for the solar-composition 
atmosphere. Nevertheless, the neutron star flux is consistent 
between the spectral models of both a blackbody and neut¬ 
ron star atmospheres. All models result in a similar goodness 
of fit. 

After the peak, at t ~ 1400 s the light curve exhib¬ 
its some variability on time-scales of ~ 10^ s. Gompared 
to the cooling model fit to the data at earlier times, the 
flux variations are both above and below the model, and 
may have a similar origin as the so-called “achromatic” 
variability, which is typically observed in some Eddington- 
limited bursts including the 1999 superburst from 4U 1820- 
30 (in’t Zand et al. 2011). In our spectral fits, the variabil¬ 
ity is strongest in the neutron star component. Since we con¬ 
strained the persistent component to be the same through¬ 
out the first orbit (Section 2.2), however, it is possible that 
the variability is actually part of the persistent emission. In¬ 
stead of the neutron star envelope, it may be more likely 
that the physical origin of the variability is found in, e.g., 
instabilities in the accretion disc. 

A systematic uncertainty in the conversion of flux to 
luminosity is introduced by anisotropy of the emission. As¬ 
suming an inclination angle of 60° (e.g., Pandel et al. 2008) 
and a thin fiat disc, the neutron star’s intrinsic luminosity is 
1/3 larger than the value we obtained under the assumption 
of isotropic emission (Fujimoto 1988). Its main effect is an 
equivalent increase in Eig. 
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The most robust part of the superburst light curve is, 
therefore, the rise and peak, until the onset of variability. 
We, therefore, fit cooling models exclusively to this part. 

5.4 Interpretation of the spectra in the tail 

In the tail, the observed signal is weaker, and it is challenging 
to separate and interpret the spectral components. Espe¬ 
cially the spectral fits with atmosphere models are problem¬ 
atic, as at certain times the fits cannot distinguish between 
very low values of T/I/Edd and more reasonable values (Fig¬ 
ure 1). Even when disregarding these unphysical solutions, 
for all spectral models the neutron star envelope’s temperat¬ 
ure appears to drop faster than predicted by cooling models 
(Figure 2). Schatz et al. (2014) pointed out that an URCA 
neutrino cooling process could operate in the outer crust, 
depending on the mixture of nuclei present, and suggested 
that the cooling could affect the shape of superburst tails. 
However, the drop in temperature that we observe here is 
much faster than expected even with URCA cooling. The 
disagreement between the model and data occurs immedi¬ 
ately after the peak, implying that the URCA cooling source 
would have to be at a very shallow depth close to the car¬ 
bon ignition depth, which is unlikely. Furthermore, fully self- 
consistent Id models of superbursts have similar problems in 
simultaneously fitting the rise and tail (Keek & Heger 2011, 
albeit this comparison was against a lightcurve that did not 
account for increased persistent emission nor reflection), in¬ 
dicating that the mismatch is not due to, e.g., late-time re¬ 
sidual burning. The strong temperature drop is found for 
the different constraints on the normalization of the neut¬ 
ron star spectral component (BB and BB^ac) as well as 
for the combined neutron star and reflection components 
(Keek et al. 2014a). If the neutron star’s emitting area is 
kept constant, this directly translates into a fast decline of 
Fns (BB), whereas with a constant reflection fraction Fns 
more closely matches the best-fitting cooling model (BBfrac; 
Figure 4). 

The slower decay of Ens for BBfrac requires a radius 
increase of a factor 2.27 ± 0.04 in the tail, even though ra¬ 
dius expansion is typically expected to occur only at larger 
luminosities near the Eddington limit. Alternatively, devi¬ 
ations of the neutron star spectrum from a blackbody have 
manifested itself as increased normalizations in the tail of 
short X-ray bursts (e.g., Suleimanov et al. 2011). By fitting 
with atmosphere models (SI and SOOl), we confirm that this 
plays no role in our case, because we fix the normalization 
in the tail. Furthermore, the neutron star spectral compon¬ 
ent is consistent with a single-temperature blackbody, in¬ 
dicating that the fast temperature decrease is intrinsic to 
the neutron star. In this case the radius increase must be 
intrinsic to the star as well, and not an anisotropy effect 
(Lapidus & Sunyaev 1985; Fujimoto 1988). 

Finally, as the flux in the tail is dominated by the per¬ 
sistent component (Figure 1; see also Keek et al. 2014a, b), 
mischaracterisation of this component could influence the 
results of the analysis in this part of the superburst. For 
example, a boundary or spreading layer is thought to be 
present between the accretion disc and the neutron star at¬ 
mosphere (Inogamov & Sunyaev 1999). The fraction of the 
neutron star covered by this layer may evolve during the su¬ 
perburst. It is, however, challenging to decompose even high- 


quality spectra of persistent emission (Revnivtsev et al. 
2013). Better constraints on the persistent spectrum may 
provided by a future superburst observation with energy cov¬ 
erage extending below 1 keV, for example with XMM New¬ 
ton (Jansen et al. 2001) or NICER (Gendreau et al. 2012). 


6 CONCLUSIONS 

The 2001 RXTE/PCA superburst from 4U 1636-536 is the 
only high-quality observation of the intrinsic rise and peak of 
a superburst. We separate the direct thermal emission from 
the neutron star from evolving reflection and persistent con¬ 
tributions to the spectrum, and evaluate the robustness of 
the derived light curve with respect to the choice of spec¬ 
tral model, including neutron star atmosphere models. In 
the tail of the burst, where the count rate is lowest, the 
interpretation of the spectrum is ambiguous. The rise and 
peak, however, provide a robust net superburst light curve. 
We show how the steepness of the rise is determined by the 
slope of the temperature profile in the neutron star envelope. 
The light curve is fit with new cooling models that take this 
feature into account. The temperature slope with pressure 
is constrained to be T oc P“ with a « 0.25. This is incon¬ 
sistent with complete burning of the carbon to iron group 
elements, and implies significant incomplete burning. The¬ 
oretical calculations of the propagation of carbon flames in 
neutron star envelopes are needed to explain the observed 
slope. Finally, the observational upper limit to the recur¬ 
rence time of this superburst is 1.75 yr, but we find that in 
that time enough fuel was accumulated for two superbursts. 
This implies that the superburst directly prior to the one in 
2001 had escaped detection. 
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